On the distribution of transmission eigenvalues in disordered wires 
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We solve the Dorokhov-Mello-Pereyra-Kumar equation which describes the evolution of an ensam- 
ble of disordered wires of increasing length in the three cases (3 = 1, 2, 4. The solution is obtained 
by mapping the problem in that of a suitable Calogero-Sutherland model. In the (3 = 2 case our 
solution is in complete agreement with that recently found by Beenakker and Rejaei. 
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During last years an increasing interest has been at- 
tracted by the physics of quantum electronic transport 
in mesoscopic systems H. One of the main reasons 
for this interest lies in the high degree of universal- 
ity of some experimental observations. In particular it 
is well known that universal conductance fluctuations 
(UCF) occur in small metallic samples at low temper- 
atures. This universality suggests that UCF could be 
described by some relatively simple Hamiltonian, inde- 
pendent of the particular model or disorder realization 
and led to the construction of a Random Matrix The- 
ory (RMT) of quantum transport ^,||, in analogy to 
the Wigner-Dyson RMT for nuclear energy levels. The 
central object in this approach is the N x N transmis- 
sion matrix t of the conductor (where N is the number 
of scattering channels at the Fermi level). All the in- 
teresting physical properties can be evaluated once the 
transmission eigenvalues T% , T% • ■ ■ TV of the product W 
are known. For instance the conductance G is given by 
G = G Tr(ttt) = G £„T„ (with G = 2e 2 /h). A 
suitable parametrization of the eigenvalues Ti is given in 
terms of the Transfer Matrix M, whose eigenvalues Xi 
are related to the T t by A* = (1 - T i )jT l . The goal is 
then to find the probability distribution P({Xi}) for the 
eigenvalues A^. This is in general a difficult problem, but 
it simplifies in the quasi-one-dimensional regime (namely 
when the length L of the conductor is much larger than 
its width W). In this case it can be shown that eigen- 
value distribution P({Xi}) as a function of the length 
L (measured in units of the mean free path I: s = L/l) 
must obey the Dorokhov-Mello-Pereyra-Kumar (DMPK) 
equation 



JV 



7 



dXi 



dXt 



(1) 



where 7 = (3N + 2 — (3 and (3 6 {1, 2, 4} is the symme- 
try index of the ensemble of scattering matrices. (3 — 1 
if there is time-reversal symmetry; /3 = 2 if the time- 
reversal symmetry is broken and (3 — 4 if time-reversal is 
conserved, but spin-orbit interactions are important and 
the spin-rotation symmetry is broken. J{{X n }) denotes 
the Jacobian from the matrix to the eigenvalue space: 



J({K}) = Y[\X 3 -X l f. 



(2) 



i<3 



Eq. (Q) was derived by Dorokhov || from microscopic 
considerations (for (3 — 2) and by Mello, Pereyra, and 
Kumar || (for /3 = 1) by looking at the infinitesimal 
transfer matrix describing the addition of a thin slice to 
the wire. Their results were then generalized to (3 = 2, 4 
in Refs. §§. 

At a first glance, solving eq.(|]) seems to be an impossi- 
ble task (except, obviously, the case in which N = 1) due 
to the interaction among eigenvalues. However, recently, 
in two remarkable papers Beenakker and Rejaei suc- 
ceeded in solving it, in the particular case (3 = 2, by using 
Sutherland's transformation and mapping eq. (Q) into the 
Schrodinger equation (in imaginary time) of a set on N 
non-interacting fermions. 

Let us briefly review their result, which will be useful 
in our following discussion. The first step suggested in |sj 
is to choose a new set of variables {x n }, related to the 
eigenvalues {A„} by: X n — sinh 2 x n , which implies T n — 
l/cosh 2 x n . Then, by making the following substitution 



(3) 



with: 



£(2) = Y\ I sinh 2 Xj - sinh 2 I ^ TT I sinh 2x i\^ ■ ( 4 ) 
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eq.(^) can be mapped into a Schrodinger equation (in 
imaginary time) for the field ^ , 



%- = {H- U)9, 
as 
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(3(P - 2) ^ sinh2 2x i + sinh 2 2x, 
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(cosh2:cj — cosh2iEi) 2 
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2 , N(N-l)^-N(N-l)(N-2)^. (5c) 



By choosing (3 — 2, the remaining interaction terms 
among the x n disappear, the equation can be decoupled 
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and, by imposing ballistic initial conditions, one ends up 
with the solution: 

P({x n }, s) — C(s) J^[(sinh 2 x 3 - sinh 2 xi) TT(sinh2a;j) x 
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dke~^ tanh(^)fc 2m_1 Pi (ifc _ 1) (cosh2a;, i ) 
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where Pv{z) are the Legendre functions of the first kind. 
C(s) is a function which depends only on s. It can be 
fixed by imposing the requirement that P({x n }, s) is nor- 
malized to unity, and does not play any role. For this rea- 
son we shall neglect in the following similar multiplicative 
functions when writing P({x n }, s). 

The aim of this letter is to show that the Hamiltonian 
Ti defined above exactly coincides with that of a partic- 
ular Calogero- SutherlandJCS) model § (for a compre- 
hensive review see ref. [|10|) with sinh-type interactions 
and a potential of type Cat (see below). This model de- 
scribes N particles on a line, identified by their coordi- 
nates {xi\, i = 1 • • • N, interacting with the potential 
1/ sinh 2 (2^ — Xj). Its most relevant feature is that (un- 
der particular conditions discussed below, see eq.(|l^)) it 
has N commuting integrals of motions, it is completely 
integrable and its Hamiltonian can be mapped into the 
radial part of a Laplace-Beltrami operator on a suitable 
symmetric space. We shall see below that in our case 
this Laplace-Beltrami operator is related in a simple way 
to the DMPK operator D of eq. ([[]) . This allows to write 
the eigenvalue distribution P({x„}, s) as a superposition 
"zonal spherical functions" of the symmetric space. By 
using the asymptotic expansion of these functions one can 
eventually obtain the explicit expression for the eigen- 
value distribution both in the metallic and in the insu- 
lating regimes. 

In the original formulation of the CS model, the inter- 
action among the particles was simply pairwise J^]. But 
it was later realized that the complete integrability of the 
model had a deep group theoretical explanation, that the 
simple pairwise interaction was the signature of an un- 
derlying structure: namely the root lattice of the Lie al- 
gebras Am and that all the relevant properties (complete 
integrability, mapping to a Laplace-Beltrami operator of 
a suitable symmetric space) still hold for potentials con- 
structed by means of any root lattice canonically associ- 
ated to a simple Lie algebra |l0| . 

Let us call V the N dimensional space defined by the 
coordinates {xt} and x = (x%, ■ ■ -Xn) a vector in V. Let 
R = {a} be a root system in V, and R + the subsys- 
tem of positive roots of R. Let us denote with x a the 
scalar product (a, x). Then the general form of the CS 
Hamiltonian is 



H = — 
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sinh (x a ) 
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where the couplings g 2 are the same for equivalent roots, 
namely for those roots which are connected with each 
other by transformations of the Coxeter group W of the 
root system. 

The CS models which are relevant for our problem are 
those constructed by means of the Cn lattice. In this 
case the root system is R = {±2e.;, ± ± ej, i ^ j}, 
(where {e^, • • • ejv} denote a canonical basis in the space 
R n ) and the Coxeter group W coincides with the product 
of the permutation group and the group of transforma- 
tions which change the sign of the vectors {e^}. The 
corresponding Hamiltonian is: 
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sinh (xi — Xj) sinh (xi + xj) 



(8) 



By using simple identities among hyperbolic functions 
eq.(||) can be rewritten as follows: 
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T* 1 sinh (2xi 
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sinh 2 (2xi) + sinh 2 (2xj) 
(cosh(2xi) — cosh(2xj)) 2 



(9) 



with c an irrelevant constant. By setting g 2 = —1/2 and 
g\ = f3{(3 — 2)/4 we see that eq. (|8]) coincides (apart from 
the overall factor 1/7) with Ti in eq.(pb|). 

The CS hamiltonian (^J) can be mapped into the radial 
part of a Laplace-Beltrami operator B (see for instance 
Appendix D of ref. EQ]) of a suitable symmetric space 



H = £(x) 



with 
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(11) 
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and p a constant term which we shall neglect in the fol- 
lowing. At this point it is easy to relate the DMPK 
operator of eq. (0) to the B operator: 



D = ±{ax)} 2 B{t;(xT 2 , 
2 7 



(12) 



and reobtain in this way the DMPK equation. This map- 
ping was discussed in great detail in ref. ||. Notice that, 
as it was stressed in ||, it is not obvious that Sutherland's 
mapping should work also for the non-translationally in- 
variant interaction of eq.(^). It is exactly the Cn struc- 
ture underlying eq. (^]) which allows such a map flic| . 

All the irreducible symmetric spaces of classical type 
can be classified with essentially the same techniques 
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used for the Lie algebras. They fall into 11 classes la- 
belled by the type of root system and by the multiplicities 
of the various roots jll| . 

A key role in the identification of the DMPK equation 
as the radial part of a Laplace-Bcltrami operator on a 
symmetric space is played by the constants g\. In fact it 
can be shown jl0| that such an identification is possible 
only if 



2 m a (m a -2) 2 

g a = z M 



(13) 



where \a\ is the length of the root a and m a its multi- 
plicity. In our case we have m a — (3 for the short roots 
(those of the type {±ei ± e^}) and m a = 1 for the long 
roots (those of the type {2ei}). Remarkably enough these 
values coincide with the multiplicities of irreducible sym- 
metric spaces exactly for those values of (3 6 {1,2,4} 
which are physically relevant. The identification is as 
follows: 



(3 = 1 
(3 = 2 
f3 = 4 



Sp(N,K)/U(N) 
SU(N,N)/S(U(N) ® U(N)) 
SO*(AN)/U(2N) 



It is interesting to notice that this same identification 
was obtained by Huffmann |l^] by directly looking to the 
symmetry properties of the DMPK equations for various 
values of (3. 

According to eq. (|l2|) if <&k{x), x — {x\, ■ ■ ■ , xjv}, k — 
{k%, ■ ■ ■ , kjsr} is an eigenfunction of B with eigenvalue fc 2 , 
then £(x) 2 &k(x) will be an eigenfunction of the DMPK 
operator with eigenvalue k 2 /(2j). These eigenfunctions 
of the B operator are known in the literature as "zonal 
spherical functions" . In the following we shall use three 
important properties of these functions (see jl3| ) . 

1] By means of the zonal spherical functions one can 
define the analog of the Fourier transform on symmetric 
spaces: 



f(x)= / 



dk 



c{k)\ 
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(where we have neglected an irrelevant multiplicative 
constant) and in the three cases which are of interest 
for us: 



\c{k)\ 
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with 
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where T denotes the Euler gamma function. 



2] for large values of x, $fc(x) has the following asymp- 
totic behaviour: 



Wife* 1 *™) ' 
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where rk is the vector obtained acting with r S W on k. 
The important feature of eq. ( [l7| ) is that it is valid for all 
values of k. 

3] in the case (3 = 2 the explicit form of $fe(x) is 
known 0l|: 



det [Q{ 



where the matrix elements of Q are: 

(\ 1 2 

and F(a, 6, c; z) is the hypergeometric function. 
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Eq.s( J12| , [l4f|l6| ) allow to write the s-evolution of 
P({x n }, s) from given initial conditions (described by the 
function fo(k)) as follows: 



P({x n }, S ) = {ax)} 2 / f (k)e-- s <S> k (x) 



dk 



\c(k)\ 



(20) 



By inserting the explicit expression of |c(fc)| 2 and by us- 
ing the identity: 



r(zf) 



k ixk 
— — tanh — — 
2 2 



(21) 



we end up with the following general expression for 
P({i„}, s) with ballistic initial conditions (which, due to 
the normalization of ^^(x), simply amount to choosing 
fo(k) — const): 

P({x n },s) = Ux)f J dfce-^ s ^^n^ tanh (^) 

(22) 

This is the main result of this letter. Let us try now a 
few applications. 

First, as a consistency check of the whole approach, 
if we insert, in the (3 = 2 case, the explicit expression 
for $fe(a;), given in eq.s( ^8|JlS| ), into eq.(f22|) and use the 
identity 



P„(«)=J , (-i/,i/+l,l;(l-2)/2) 



(23) 



we exactly obtain (as expected) the solution, eq. (g) found 
by Beenakker and Rejaei. 

Second, if x is large (and in our framework this means 
x 2 > (2s) /j) we may insert the asymptotic expansion 
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into eg. (p2|) . The resulting behaviour of P({x n },s) 
will depend on the chosen (metallic or insulating) regime 
for k. In both cases, by setting j3 = 2 in eq.s ( |25| ) and 
( p8| ) below, we exactly recover the results described in || , 
in these regimes. Let us look at the two cases separately. 

Insulating regime (k <C 1). 
In the k — > limit the T functions in eq.(jl^) can be 
approximated according to: 



+ Q iy , f3 e {1,2,4} 



(24) 



Then, by rewriting both the product ~[ i<: j(k 2 — kj), and 
the sum over the exponentials in (Jl7) as determinants, 
the integration over k becomes straightforward and gives: 



P({x n }, s) — |sinh 2 Xj — sinh 2 Xi\ 2 [( 



2 _ 2 



)] 



x ]J |^exp(-a; 2 7/(2s))a; 4 (sinh2a; 4 ) 1/2 



(25) 



Ordering the x„'s from small to large and using the 
fact that in this regime 1 <C x\ <C x% <C • • • "C xn we 
can approximate the eigenvalue distribution as follows: 



P({x n },s) = Y[exp[-(j/(2 S ))(i 



(26) 



where x n = 1)), in agreement with the result 

obtained by Pichard, |l5| by directly solving the DMPK 
equation in this regime. 

Metallic regime (k ^> 1). 
In this case one must use the asymptotic expansion: 



The integration over fc is less simple in this case and, (to 
be consistent with the regime of validity of eq.(|l7|)) in the 
resulting expression only the highest powers of {x^f^-) 
must be taken into account. We find: 



values of k (see sect. 8 of ref. ju| and ref. §J). The 
coefficients of this series can be constructed recursively, 
thus allowing to study the behaviour of eq. (^2j) even for 
values of x smaller than the above mentioned threshold. 
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P({x n }, s) = J^J |sinh 2 Xj — sinh 



2 I T 12 



„2|l 



i<j 



Y[ [exp(-x 2 7 /(2 S ))(x l sinh2x l ) 1 / 2 ] . (28) 



In agreement with the exact result of |g] for (3 = 2 and 
with the /3 dependence found by Chalker and Macedo |lr| 
through a direct integration of the DMPK equation. 

According to ref. [ flO[ and ref. |l~3| , the regime of va- 
lidity of eq.s(p5|) and (|28|) is x 2 > (2s) /j, which gives in 
in the large N limit x 2 > (2s)/ (/37V). Notice however, as 
a concluding remark, that eq.(|l7|) is only the first term 
of a series which converges absolutely to $fe(x) for all 



4 



